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1 . Introduction 

This project has two goals: First, to develop very fast algorithms 

for solution of the transient and stationary-state equations describing 
combustion chemical kinetics in batch and stirred-tank chemical reactors, 
for imbedding in larger hydrodynamic codes which model mixing, heat 
transfer and chemical change in gas turbine combustors and other combus- 
tion devices. The second goal is to develop "user-friendly" interactive 
computer programs which combine the fast kinetics algorithms with graphics 
and user selection menus to create a dynamic environment for parametric 
variation, preliminary design, and in general permitting the user to gain 
an improved understanding of the thermal -chemical -phys ical behavior of 
combustion processes, without being hampered by considerations of numerical 
methods, stability, convergence or computational efficiency. 

2. Batch Kinetics (l-D) Algorithm Development 

A typical batch combustion problem at constant pressure consists of 
three distinctly different chemical-physical regimes, each of which 
requires a unique algorithm for the peculiar form of stiffness of the 
governing set of ordinary differential equations. (The term "stiff", 
coined by Hi rschfelder , means that the characteristic time constants of 
the coupled set of equations differ by many orders of magnitude, so that 
integration of the equations by conventional methods would require time- 
steps very much smaller than desired, resulting in excessively large 
computational work. [1]) 

A single-step integration algorithm was developed which automatically 
identifies equations as having either positive or negative time constants, 
and then approximates the unstable equations (positive time constants) 
with the trapezoidal rule, and the stable equations (negative time 
constants) with a decaying exponential function. This "exponent ial -fi tted 
trapezoidal rule" was adapted form the work of Liniger and Willoughby [2] 
and of Brandon [3] • 

During induction and early heat-release regime, the species equations 
are dominated by positive time constants, and the temperature also 
exhibits a positive time constant. Since very smal 1 steps are required 
for integration of unstable equations, a simple predictor-corrector scheme 
with functional iteration assures the least computational work possible. 
However, during late heat release and equilibration, when the temperature 
and species equations exhibit negative time constants, large stepsizes 
are now possible, so that Newton-Raphson iteration with calculation of 
the full Jacobian matrix is the optimal convergence method. At the present 
level of development, the batch kinetics code CREK-ID (for "combustion 
reaction kinetics, one-dimensional") appears to be at least five times 
faster than LSODE, the Lawrence Livermore Laboratory implementation of 
the Gear-B algorithm [4]. Further details of the algorithm are given 
in Reference 5. 
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3. Stirred Reactor (O-D) Algorithm Development 


An earlier packaged code, CREK [6], solved for time-stationary 
states of a perfectly-stirred chemical reactor. A major improvement 
has been made in the automatic estimate-generating feature of the CREK 
code. In this approach, the transient sti rred-reactor equations are 
solved for the predictor step only, using the same techniques developed 
for the batch kinetics CREK- ID code. The resulting approximate solution 
gives a very good estimate of the steady-state solution, which is then 
solved for exactly using the original CREK algorithm. As a result, the 
speed and reliability of the "CREK-OD" program appears to have been 
increased by at least an order of magnitude over the original CREK code. 

4. Interactive Computer Programs 

Two "user-friendly”, menu-driven programs have been written to 
date. The first program, EQLBRM, solves for chemical equilibrium states 
for combustion of arbitrary hydrocarbon fuels with air, at a specified 
constant pressure, fuel-air ratio and initial temperature of fuel and 
air. The second program, AVC0-M2, is a reactor-theoretic simulation of 
a gas turbine combustor [7]. The combustor is assumed to consist of 
up to nine "flow elements" connected in series, and optionally one flow 
element connected in recycle. Each flow element may be a zero-dimensional 
or perfectly-stirred reactor (PSR) , a one-dimensional or plug-flow 
reactor (PFR) , or a non-reacting mixer (MIX). The user is prompted to 
select the order, type and size of the flow elements, the recycle path, 
and the distribution of fuel and air into each individual flow element. 

On demand, the program also displays a crude graphic schematic of the 
combustor model under consideration. 
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Variation with time of temperature and chemical species 
mole fraction in adiabatic batch reaction. 
(Pyrolized CH^-air, Stoichiometric, T=1000K, P=10 atm) 
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Conventional staggered finite-difference grid for calculation 
within a flow domain. Points H (“high") and L ('-'low") are 
above and below the plane of E-W-N-S-P. PP refers to condi- 
tions at P at the previous time step. 
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HOMOGENEOUS, ADIABATIC, GAS-PHASE 
REACTION AT CONSTANT PRESSURE 
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Batch kinetics (1-D) 
Stirred-reactor kinetics (O-D) 
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"LOCALLY EXACT" SOLUTIONS 
OF THE TRANSIENT EQUATIONS 


1 . Decoupl i ng 




A. . B.a. 


2. Integration 


n+1 


“i 




a" . 

“1 ^ B7 




"+1 

or a, ' Of + hfj 


B,h 


where h = t'^''’^ - t^ 


3. Analysis 

The locally exact solution is stable (for h — > “) if and 
only if B^.< 0. 
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FORMALLY DECOUPLING THE 
EOUATIONS: "L - FORMULATION" . 


let 






- Q. 




then 


dt 


q, - L,a. . i (o, - a,) 


Q. tio* 
pV 1 
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''i = 
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ni * 

Q • “*■ “TT O . 
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. . .A^ is positive-definite; negative-definite 
*. "L-formulated" equations are always 
stable , but may be inaccurate. 
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FUNCTIONALLY DECOUPLING THE 
EQUATIONS : "Z-FORMULATION" 


let 


da. 


= f.. = 


dt 'i 


xO . 7O / 0\ 

f i * Z, (o, - a, ) 


where Z. e 


df. 


i da^ 


Z + -?r 


k=1 H 




(and where t = - Z h|^f|^ / Z C„ a^) 


k Pk 


da. 


then 


' - [f? - Z°a°] + Z°a. 


dt 


=> A. = 


[f“ 


7O On 

Z,o,] 


6 . 

1 




. . . neither A^. nor are positive-definite or negative-definite. 

"Z-formulated" equations are always accurate , but may be unstable . 
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HARK2I )K>Jo»c 

■MARK2I* IS AN INTERACTIVE VERSION OF THE MARK-II COMBUSTOR MODEL* 
THIS IS A PRELIMINARY DESIGN TOOL, A MEANS OF GAINING INTUITIVE INSIGHT 
INTO EFFECTS OF CHANGES IN FUEL-AIR MIXING OR PARTITIONING ON TURN-DOWN 
RATIO, COMBUSTION EFFICIENCY AND POLLUTANT FORMATION RATES* 

AN INITIAL DATA SET IS TAKEN FROM DATA FILE -MARK2*DAT- BUT CAN BE ALTERED 
INTERACTIVELY, AND USED IN CONSECUTIVE RUNS* 

MAF^K-II REPRESENTS A SIMPLE BRAGG COMBUSTOR CONSISTING OF A 
MAXIMUM OF 9 FLOW ELEMENTS WITH THE ADDITION OF A SINGLE RECYCLE 
ELEMENT* FLOW ELEMENT TYPES MAY INCLUDE: 

1) NON-REACTING MIXERS CMIX*), IN WHICH THE CHEMICAL REACTIONS ARE 
ASSUMED TO HAVE STOPPED DURING THE MIXING PROCESS; 

2) PERFECTLY STIRRED REACTORS CPSR*), WITHIN WHICH INTENSE SELF- OR 
DACK-MIXING IS ASSUMED TO OCCUR, SO THAT THERE ARE NO AXIAL GRADIENTS; 

3) PLUG FLOW REACTORS CPFR*)* 

THE USER MAY DEFINE THE MODEL AS HAVING UP TO 9 ELEMENTS IN 
SERIES WITH AIR AND FUEL INLET JETS AT EACH ELEMENT* THE RECYCLE ELEMENT 
MAY DE OF ANY OF THE THREE FLOW TYPES, AND MUST RECYCLE FROM A HIGHER 
NUMBERED ELEMENT TO A LOWER* COOLING BOUNDARY LAYER EFFECTS AND CHEMICAL 
REACTIONS WITHIN THE BOUNDARY LAYER ARE NOT CONSIDERED* 

PLEASE WAIT A MOMENT WHILE INITIALIZATION IS COMPLETED. 

— INITIALIZED ~ PRESS <RETURN> TO BEGIN — 


*** INPUT DATA 


FLOW 

AREA 

LENGTH 

FLOW 

ELEMENT* 

(SQ*IN) 

(INCHES) 

TYPE 

1 

1*4A00ET02 

4*0000E-01 

PSR 

n 

l*4600E+02 

1*0000E-01 

MIX 

3 

1*4A00ET02 

l*5000E+00 

PSR 

4 

1*4S00EF00 

2*0000E-01 

MIX 


INLET AIR 
<LBM/S) 
l*3450E+00 
6*9600E-"01 
0*0000E-01 
7*9500E-01 


INLET FUEL 
(LBM/S) 
9.7200E-0 
O.OOOOE-0 
0*0000E-01 
0*OOOOE-Ol 


AIR TEMP « 2»1000ET02 F 

FUEL TEMP 8*0000E+01 F 


COMBUSTOR PRESSURE = 2*2100E+00 ATM 

LOWER HEATING VALUE ^ 1*8500E+04 BTU/LBM 


SELECT AN OPTION BY NUMBER: 

-0- RUN WITH THIS DATA SET 
-1- CHANGE AIR TEMPERATURE 
-2“ CHANGE FUEL TEMPERATURE 
•3- CHANGE LOWER HEATING VALUE 
OPTION? (0-7) 7 


-4- CHANGE NOMINAL COMBUSTOR PRESSURE 
-5- CHANGE RECYCLE ELEMENT STATUS 
-6- CHANGE FLOW ELEMENTS STATUS 
-7- INSPECT SCHEMATIC MODEL LAYOUT 


MARK-II MODEL SCHEMATIC LAYOUT 


A F A A 

1 ! / 1 ! 

^ 

* 1 t t 2 t t Z t ^ A t 
->*PSR)|t— >*MIX* — >»PSR3|c— »»cMIX5»c — > 


PRESS <RETURN> TO CONTINUE 
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%%% INPUT DATA 


FLOW 

AREA 


LENGTH 

FLOW 

INLET AIR 

INLET FUEL 

ELEMENT* 

(SQ*IN) 


(INCHES) 

TYPE 

(LBM/5) 

(LBM/S) 

1 

1 ♦>>600ET02 


>^.0000E-01 

PSR 

l*3A50E+00 

9t7200E-02 

n 

1.4600E+02 


1*OOOOE-Ol 

NIX 

6*?600E-"01 

O.OOOOE -o:i 

3 

1.4600E+02 


l»5000E+00 

PSR 

0,0000E“01 

0*0000E-01 

A 

l*4500E+00 


2*0000E‘ 01 

MIX 

7.9S00E--01 

O.OOOOE-01 

AIR TEMP = 

2*1000ET02 

F 

COMBUSTOR 

PRESSURE 

^ 242:l00E+00 

ATM 

UEL TEMP 

OiOOOOE+Ol 

F 

LOWER HEATING VALUE 

= l*8500E+0^ 

BTU/LBM 


SELECT AN OPTION BY NUMBER; 

-0- RUN WITH THIS BATA SET 
-1- CHANGE AIR TEMPERATURE 
-2- CHANGE FUEL TEMPERATURE 
-3- CHANGE LOWER HEATING VALUE 
OPTION? <0-7> 5 


-A- CHANGE NOMINAL COMBUSTOR PRESSURE 
-5- CHANGE RECYCLE ELEMENT STATUS 
-A- CHANGE FLOW ELEMENTS STATUS 
-7- INSPECT SCHEMATIC MODEL LAYOUT 


RECYCLE ELEMENT STATUS: 

CROSS-SECTIONAL AREA - I.4600EL02 SO. INCH 

ELEMENT LENGTH ^ l.OOOOE+00 INCHES 

FLOW TYPE MIX 

RECYCLES 0.00% OF ELEMENT * 1 OUTFLOW TO ELEMENT ♦ 1 INFLOW 

SELECT AN OPTION BY NUMBER; 

0- NO MORE CHANGES TO RECYCLE ELEMENT 

1- CHANGE AREA 

2- CHANGE LENGTH 

3- CHANGE FLOW TYPE 

A- CHANGE RECYCLE INTAKE PERCENTAGE 
5- CHANGE RECYCLE PATH 
ELIMINATE RECYCLE 

OPTION? <0-6) A 

ENTER PERCENT OF OUTFLOW TO BE RECYCLED: 20, 


RECYCLE ELEMENT STATUS; 

CROSS-SECTIONAL AREA “ 1.4600E+02 SO. INCH 

ELEMENT LENGTH - l.OOOOE+00 INCHES 

FLOW TYPE - MIX 

RECYCLES 20.00% OF ELEMENT * 1 OUTFLOW TO EL.EMENT ♦ 1 INFLOW 

SELECT AN OPTION BY NUMBER; 

0- NO MORE CHANGES TO RECYCLE ELEMENT 

1- CHANGE AREA 

2- CHANGE LENGTH 

3- CHANGE FLOW TYPE 

4- CHANGE RECYCLE INTAKE PERCENTAGE 

5- CHANGE RECYCLE PATH 

6- ELIMINATE RECYCLE 

OPTION? (0-6) 5 

ENTER ELEMENT WHOSE OUTFLOW IS TO BE RECYCLED: (1- 4) 3 
ENTER ELEMENT WHOSE INFLOW IS TO RECEIVE RECYCLE; (1-3) 2 
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RECYCLE ELEMENT STATUS: 

CROSS-SECTIONAL AREA 1.^600E+02 SQ/INCH 
ELEMENT LENGTH = l,OOO0E+OO INCHES 

FLOW TYRE == MIX 

RECYCLES 20*00% OF ELEMENT ♦ 3 OUTFLOW TO ELEMENT ♦ 2 INFLOW 

SELECT AN OPTION BY NUMBER: 

0- NO MORE CHANGES TO RECYCLE ELEMENT 

1- CHANGE AREA 

2- CHANGE LENGTH 

3- CHANGE FLOW TYPE 

4- CHANGE RECYCLE INTAKE PERCENTAGE 

5- CHANGE RECYCLE PATH 

6- ELIMINATE RECYCLE 

OPTION? (0-6) 0 


INPUT 

DATA *** 






FLOW 

AREA 

LENGTH 

FLOW 


INLET AIR 

INLET FUEL 

ELEMENTI 

(SQ*IN) 

(INCHES) 

TYPE 


(LBM/S) 

(LBM/S) 

1 

l*4600E+02 

4.0000E-01 

P9R 


l*3450Ef00 

9*7200E-02 

2 

1*4600ET02 

l.OOOOE-01 

MIX 


6*9600E-01 

0*0000E-01 

3 

l*4600E+02 

1 ^SOOOETOO 

PSR 


0*OOOOE -01 

0*0000E-01 

4 

l*4500E+00 

2*0000E-0I 

MIX 


7*9500E-01 

0*0OOOE“Ol 

RECYCLE 

1*4600ET02 

l.OOOOE+OO 

MIX 


RECYCLE 2' 

0*00% OF * 3 






OUTFLOW TO 

♦ 2 INFLOW 

AIR TEMP 

= 2*1000EF02 F 

COMBUSTOR PRESSURE 

2*2100E+00 

ATM 

FUEL TEMP 

- 8*0000E+01 F 

LOWER HEATING VALUE 

l*8500E+04 

BTU/LBM 

SELECT AN 

OPTION BY number: 






•0- RUN WITH THIS DATA SET 


CHANGE 

NOMINAL COMBUSTOR 

PRESSURE 

-1- CHANGE 

AIR TEMPERATURE 

-5- 

CHANGE 

RECYCLE ELEMENT STATUS 

-2- CHANGE 

FUEL TEMPERATURE 

-6- 

CHANGE 

FLOW 

ELEMENTS STATUS 

-3- CHANGE 

LOWER HEATING VALUE -7- 

INSPECT 

SCHEMATIC MODEL 

LAYOUT 


OPTION? (0-7) 7 


MARK-II MODEL SCHEMATIC LAYOUT 


A F 
1 1 


A 

1 



A 

1 

^ ^ 






* 1 * 


* 2 * 

^ Z ^ 


* 4 

5KPSR*- 

— > 

•*MIX>f^-— 

»KPSR)1^- 


^^*:MIX* 


1 



1 

^ ^ 




PRESS <RETURN> TO CONTINUE 
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*** INPUT DATA 

FLOW AREA LENGTH FLOW INLET AIR INLET FUEL 

ELEMENT* <SQ*IN) (INCHES) TYPE (LBM/S) (LBM/S) 

1 l*4600E+02 4*0000E“01 PSR l*3450E+00 9>7200E-02 

2 l*4600E+02 1,0000E-01 MIX 6^9A00E-01 O.OOOOE-Ol 

3 1 ♦4600E+02 .1 *5000E+00 PSR 0*0000E~01 0,0000E~01 

4 l*4500E+00 2*0000E-01 MIX 7.9500E-01 O.OOOOE~01 

AIR TEMP = 2.1000E+02 F COMBUSTOR PRESSURE == 2*2100E+00 ATM 

FUEL TEMP 8*OOOOE+Ol F LOWER HEATING VALUE - l*8500E+04 BTU/LBM 

SELECT AN OPTION BY NUMBER: 

-0- RUN WITH THIS DATA SET -4- CHANGE NOMINAL COMBUSTOR PRESSURE 

-1- CHANGE AIR TEMPERATURE CHANGE RECYCLE ELEMENT STATUS 

-2~ CHANGE FUEL TEMPERATURE -6- CHANGE FLOW ELEMENTS STATUS 

-3- CHANGE LOWER HEATING VALUE -7- INSPECT SCHEMATIC MODEL LAYOUT 

OPTION? (0-7) 6 


TO CHANGE DATA, ENTER THE CODE AND ELEMENT (IE -B2* FOR *2'S LENGTH) 

CODE -A- CODE 'B- CODE “C* CODE Tr- 

FLOW AREA LENGTH FLOW INLET AIR 

ELEMENT* (SQ*IN) (INCHES) TYPE (LBM/S) 

1 1.4600E+02 4*0000E-01 PSR 1 *3450E+00 

2 l*4600E+02 1,0000E“01 MIX 6.9600E-01 

3 1 ♦ 4600E+ 02 1 ♦ 5000ETOO PSR 0 * OOOOE-01 

4 1»4500ET02 2»0000E-01 MIX 7*9500E-01 

CODE "F- “ REMOVE ENTIRE FLOW ELEMENT 

CODE 'G* - ADD ANOTHER FLOW ELEMENT (NO NUMBER NEEDED) 

(BLANK ENTRY RETURNS TO MAIN MENU) 

ENTER CODE, ELEMENT* A4 

ENTER CROSS-SECTIONAL AREA OF ELEMENT ♦ 4 (SQUARE INCHES): 2»00E2 
ENTER CODE, ELEMENT* A4 

ENTER CROSS-SECTIONAL AREA OF ELEMENT * 4 (SQUARE INCHES): 1,45E2 
ENTER CODE, ELEMENT* 


INPUT DATA 


Fl-OW 

AREA 


LENGTH 

FLOW 

IN’LET AIR 

INLET FUE 

ELEMENT* 

(SQ* IN) 


(INCHES) 

TYPE 

(LBM/S) 

(LBM/S) 

1 

1»4600E+02 


4,0000E 01 

PSR 

1 ,34S0E+00 

9.7200E-02 

r> 

l*4600E+02 


1*0000E-01 

MIX 

6*9600£-01 

O.OOOOE-01 

3 

1*4600E+02 


1 tSOOOEiOO 

PSR 

0*0000E-01 

O.OOOOE-01 

4 

l*4500E+00 


2.0000E--01 

MIX 

7.9500E-01 

0,OOOOE-Ol 

RECYCLE 

l,4A00E+02 


1*0000ET00 

MIX 

RECYCLE 20 
OUTFLOW TO 

.00% OF * 3 
* 2 INFLOW 

AIR TEMP = 

2*1000ET02 

F 

COMBUSTOR 

PRESSURE 

2*2100E+00 

ATM 

FUEL TEMP 

8.0000E+01 

F 

LOWER HEATING VALUE == 

1*8500ET04 

BTU/LBM 


SELECT AN OPTION BY NUMBER: 

-0- RUN WITH THIS DATA SET 
-1- CHANGE AIR TEMPERATURE 
-2- CHANCE FUEL TEMPERATURE 
-3- CHANGE LOWER HEATING VALUE 
OPTION? (0-7) 2 

ENTER FUEL TEMPERATURE (DEG FAHRENHEIT): 90* 


-4- CHANGE NOMINAL COMBUSTOR PRESSURE 
-5- CHANGE RECYCLE ELEMENT STATUS 
-6- CHANGE FLOW ELEMENTS STATUS 
-7- INSPECT SCHEMATIC MODEL LAYOUT 


CODE -E* 
INLET FUEL 
(LBM/S) 
9,7200E-02 
0,OOOOE-Ol 
0*OOOOE-Ol 
0»OOOOE-O:l 
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JK** INPUT DATA 


FLOW 

AREA 

LENGTH 

FL.OW 


INLET AIR 

INLET FUEL 

ELEMENT* 

(SQ*IN) 

(INCHES) 

TYPE 


(LBM/S) 

(LBM/S) 

1 

l,4600E+02 

4,0000E--01 

PGR 


l*34fJ0E+00 

9*7200E-02 

2 

1.4600E+02 

1*0000E“01 

MIX 


6.9600E-01 

O.OOOOE'OL 

3 

1.4600E+02 

l*5000E+00 

PBR 


0*0000E“01 

0*0000E-Ol 

4 

1*4500ET00 

2*0000E-01 

MIX 


7*9500E“01 

O.OOOOE-01 

r^ECYCLE 

l*4600E+02 

l^OOOOETOO 

MIX 


RECYCLE 20.00% OF ♦ 3 






OUTFLOW TO 

♦ 2 INFLOW 

(MR TEMP =: 

2.1000E+02 F 

COMBUSTOR PRESSURE 

2*2100E+00 

ATM 

FUEL TEMP - 

9»0000E+01 F 

LOWER HEATING VALUE = 

1 .8500E+04 

BTU/LBM 

SELECT AN OPTION BY NUMBER: 






-0- RUN WITH THIS DATA SET 

.-.4- 

CHANGE 

NOMINAL COMBUSTOR 

PRESSURE 

-1- CHANGE 

AIR TEMPERATURE 

-5- 

CHANGE 

RECYC 

LE ELEMENT STATUS 

-2- CHANGE 

FUEL TEMPERATURE 

-6- 

CHANGE 

FLOW 

ELEMENTS STATUS 

-3- CHANGE 

LOWER HEATING VALUE -7- 

INSPECT 

SCHEMATIC MODEL 

LAYOUT 


OPTION? (0-7) 0 

RUN WITH THIS DATA SET ARE YOU SURE? (Y/N) Y 


t** output: *** 

TEMP AND COMPOSITION FOR MARK-II MODEL: 

PRESSURE = 2,21 ATM OVERALL E.R. = 0.4995 

FUEL FLOW = 349,92 LBM/HR AIR FLOW = 2.84 LBM/SEC 

20.0 PCT. OF FLOW RECYCLED FROM OUTLET OF 
ELEMENT NO. 3 TO INLET OF ELEMENT NO, 2 

OTHER output: SELECT BY NUMBER - 

0- OUIT /BEGIN NEXT INPUT DATA SET 

1- INSPECT MOLE FRACTIONS 

2- INSPECT MASS FRACTIONS 

3- INSPECT EMISSIONS INDEX 

4- INSPECT MISC. OTHER OUTPUT 


5- INSPECT 
OPTION? (O' 

CURRENT INPUT 
-4) 1 

BATA SET 





MOLE FRACTIONS X(I) 






SPECES 


unt 

f2» 

5»C3* 

)»C43|£ 

EQL 

RCY 

C12H23 


7*484D-05 

4*070D-05 

2*1540-07 

1*5720-07 

2*8860-19 

2*1540- 07 

C2H2 


l*080D-02 

6*32211-03 

2*3000-03 

1*6790-03 

2*8860-19 

2*3000-03 

C2H3 


9*7a2D-04 

5*3270-04 

6*5580-06 

4*7850-06 

2*8860-19 

6*5580-06 

C2H4 


7*447D“04 

4*054D“04 

4*6250-06 

3*3750-06 

2*8860-19 

4*6250-06 

CH20 


5»914D“05 

3*2980-05 

4*3150-06 

3*1490-06 

1*6960-18 

4*3150-06 

CH3 


5*570D“05 

3*0920-05 

3*3320-06 

2*4310-06 

2*8860-19 

3*3320-06 

CH4 


3*81511-12 

2*5900-12 

2*6930-12 

1*9650-12 

2*8860-19 

2*6930-12 

CO 


1 *428D-02 

8*3350-03 

2*9020-03 

2*1180-03 

1*6940-06 

2*9020-03 

C02 


9*676Ei“02 

6*9700-02 

8*6160-02 

6*2880-02 

6*6580-02 

8*6160-02 

HCO 


l*064H-04 

5*8480-05 

3*3660-06 

2*4560-06 

5*7520-15 

3*3660-06 

H 


l*246D-03 

7*2750-04 

2*5540-04 

1*8640-04 

1*8140-08 

2*5540-04 

H2 


2*892D-03 

1*6810-03 

5*5690-04 

4*0640 04 

5*7610-07 

5*5690-04 

H20 


1*108D-01 

7*7160-02 

8*5410-02 

6*2330-02 

6*5690-02 

8*5410-02 

H02 


2*86111-04 

l*724ti-04 

8*5720-05 

6*2550-05 

7*4800-08 

8.5720-05 

N 


3.167D-07 

1.7510-07 

3*4790-08 

2*5390-08 

1*6020-13 

3*4790-08 

NO 


3*253D-05 

4*9410-05 

5*8140-05 

4*2430-05 

1*2010-03 

5*8140-05 

N02 


l*374D-07 

2*8790-07 

4*4420-07 

3*2410-07 

5*2850-06 

4*4420-07 

N2 


7*249D-01 

7*4690-01 

7*5120-01 

7*6170-01 

7*6340-01 

7*5120-01 

0 


1*53111-03 

1*0220-03 

9*5940-04 

7*0010-04 

1*9120-06 

9*5940-04 

OH 


4*489D-03 

2*9650-03 

2*6510-03 

1*9350-03 

7*1520-05 

2*6510-03 

02 


3*000D-02 

8*3860-02 

6*7440-02 

1*0600-01 

1*0100-01 

6*7440-02 

PRESS <return: 

> TO CONTINUE 
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output: 


TEMP AND COMPOSITION FOR MARK-II MODEL: 

PRESSURE ^ 2f21 ATM OVERALL E.R* = 0,4995 

FUEL FLOW = 349,92 LBM/HR AIR FLOW = 2,84 LBM/SEC 

20,0 PCT, OF FLOW RECYCLED FROM OUTLET OF 
ELEMENT NO, 3 TO INLET OF ELEMENT NO. 2 

OTHER output: SELECT BY NUMBER - 

0- QUIT / BEGIN NEXT INPUT DATA SET 
1“ INSPECT MOLE FRACTIONS 

2- INSPECT MASS FRACTIONS 

3- INSPECT EMISSIONS INDEX 

4- INSPECT MISC, OTHER OUTPUT 

5- INSPECT CURRENT INPUT DATA SET 
OPTION? (0-4) 4 


ELEMENT NO, *1* *2# *3* 

ELEMENT TYPE PSR MIX PSR 

EQUIV RATIO l,053D+00 A,941D-01 6,941D-01 
RES TIMEf SEC 5.647D-04 9,558D-05 1.261D-03 
AREA, SQ IN l,460D+02 1.460D+02 l,460D+02 
VELO, FT/SEC 5.903D+01 8,719D+01 9.909D+01 
FLOW, LBM/SEC l,442D+00 2,673D+00 2,673D+00 
AXIAL LOC, IN 4,000D-01 5,000D-01 2.000D+00 
ENTH , BTU/LBM-2 ♦ 219D+01-4 , 426D+00-4 , 427D+00 
TEMP EFF 8.192-01 7.850D-01 9.317D-01 

TEMP, DEG F 3.105D+03 2,400D+03 2.810D+03 
PRESS <R£TURN> TO CONTINUE 


*4* 

MIX 

4.995D-01 
1.879D-06 
1 ,45011+00 
8.872D+03 
2.933D+00 
2.200D+00 
5.530D+00 
9.260D-01 
2,192D+03 


EQL 

EQL 

4.995D-01 

O.OOOD-01 

0,000D-01 

O.OOOD-01 

2.933D+00 

2.200D+00 

5,530D+00* 

1 .OOOD+00 

2.350D+03 


output: 

TEMP AND COMPOSITION FOR MARK-II MODEL: 

PRESSURE == 2,21 ATM OVERALL E.R. = 0.4995 

FUEL FLOW == 349,92 LBM/HR AIR FLOW = 2.84 LBM/SEC 

20.0 PCT. OF FLOW RECYCLED FROM OUTLET OF 
ELEMENT NO. 3 TO INLET OF ELEMENT NO. 2 

OTHER output: SELECT BY NUMBER - 

0- QUIT / BEGIN NEXT INPUT DATA SET 

1- INSPECT MOLE FRACTIONS 

2- INSPECT MASS FRACTIONS 

3- INSPECT EMISSIONS INDEX 

4- INSPECT MISC. OTHER OUTPUT 

5- INSPECT CURRENT INPUT DATA SET 
OPTION? <0-4) 3 


RCY 
RCY-MIX 
A.941D-01 
4.205D-03 
1.460D+02 
1 .982D+01 
5*346D-01 
4,304D+10 
4,380D+00 
9.317D-01 
2,81011+03 
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EMISSION 

INDEX EI(I) , GM 

i I/KG FUEL 




SPECES ■! 

■ *1* 


)|c3* 


EQL 

RCY 

C12H23 

6.545D+00 

6.552D+00 

3*4490-02 

2*7590-02 

5*0490-14 

6*8970-03 

C2H2 

l*470D+02 

l*5B4D+02 

5*7310+01 

4*5850+01 

7*8580-15 

1*1460+01 

C2H3 

1*38311+01 

1*386D+01 

1*6970-01 

1*3580-01 

8*1620-15 

3*3940-02 

C2H4 

1*092D+01 

1*095D+01 

1. 2410 01 

9*9300-02 

8*4660-15 

2*4830-02 

CH20 

9*282D-01 

9*530D-01 

1*2400-01 

9*9130-02 

5*3250-14 

2*4790-02 

CH3 

4*377D-01 

4*47311-01 

4.7930-02 

3*8340-02 

4*5370-15 

9*5860-03 

CH4 

3. 199D-08 

3*999D-08 

4.1340-08 

3.3070-08 

4*8410-15 

8*2670-09 

CO 

2*091D+02 

2*246D+02 

7*7780+01 

6*2220+01 

4*9630-02 

1*5560+01 

C02 

2*226D+03 

2*952D+03 

3.628Dt-03 

2*9020+03 

3*1560+03 

7*2560+02 

HCO 

l*A14D+00 

l*633D+00 

9 *3440-02 

7*4750-02 

1*7450-10 

1*8690-02 

H 

6*565D-01 

7*056D-01 

2*4630-01 

1*9710-01 

1*9120-05 

4*9260-02 

H2 

3*047D+00 

3*262D+00 

1*0740+00 

8*5920-01 

1.2140-03 

2*1480-01 

H20 

l*043D+03 

l*338D+03 

1*4720+03 

1*1780+03 

1*2380+03 

2*9440+02 

H02 

4*936D+00 

5*476D+00 

2*7070+00 

2.1660+00 

2*5820-03 

5*4140-01 

N 

2*31BD-03 

2*361D-03 

4.662D-04 

3*7300-04 

2*3460-09 

9*3250-05 

NO 

5*102D-01 

1*42711+00 

1*6690+00 

1*3350+00 

3*7700+01 

3*3380-01 

N02 

3*303D-03 

l*275D-02 

1*9550-02 

1*5640-02 

2*5420-01 

3*9100-03 

N2 

1.061D+04 

2*013D+04 

2*0130+04 

2*2380+04 

2*2360+04 

4.0270+03 

0 

1.281D+01 

1*574D+01 

1. 4690+01 

1*1750+01 

3.1980-02 

2*9370+00 

OH 

3*990D+01 

4*853D+01 

4*3140+01 

3.4520+01 

1*2720+00 

8*6290+00 

02 

5*018D+02 

2*5820+03 

2*0650+03 

3*5570+03 

3*3800+03 

4.1300+02 

PRESS <RETURN> TO CONTINUE 






output: 

TEMP AND COMPOSITION FOR MARK-II MODEL: 

PRESSURE = 2*21 ATM OVERALL E*R* = 0*4995 

FUEL FLOW = 349*92 LBM/HR AIR FLOW 2*84 LBM/SEC 


20*0 PCT* OF FLOW RECYCLED FROM OUTLET OF 
ELEMENT NO* 3 TO INLET OF ELEMENT NO* 2 


OTHER output: SELECT BY NUMBER - 

0- QUIT / BEGIN NEXT INPUT DATA SET 

1- INSPECT MOLE FRACTIONS 

2- INSPECT MASS FRACTIONS 
3“ INSPECT EMISSIONS INDEX 

4- INSPECT MISC. OTHER OUTPUT 

5- INSPECT CURRENT INPUT DATA SET 
OPTION? (0-4) 0 

no YOU WANT A HARD COPY OF THE INPUT AND OUTPUT DATA FOR THIS RUN? (Y/N) Y 


DO YOU WANT ANOTHER RUN? (Y/N) N 
QUIT ARE YOU SURE? (Y/N) Y 


— FORTRAN STOP 
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PRESENT STATUS OF PROJECT - OCTOBER 1982 


Aval lable now: 

IMPROVED CREK (-0D) 

Including CREK preprocessor routines, 
sample data dekc, 1976 CREK documentation 

EQLBRM and MARK21 

interactive programs 

CREK- ID 

with CREK preprocessor routines and 1982 paper 
FOR PRELIMINARY EVALUATION AND CRITICISM ONLY! 


ALL SOURCE CODES; ANS I -standard 1966 FORTRAN-IV, ASCII 
MEDIA: Listing, or 8” floppy disk, single density, 

IBM 37^0 formatted; DEC RXOl in RT-11 or FILES-11 
other media negotiable 
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